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Abstract A method of calculating induced electric field is presented in this paper. Induced 
electric field in solar atmosphere is derived by the time variation of magnetic field when 
the charged particle accumulation is neglected. In order to get the spatial distribution 
of magnetic field, several extrapolation methods are introduced. With observational data 
from Helioseismic and Magnetic Imager (HMI) aboard the NASA's Solar Dynamics 
Observatory (SDO) on May 20th, 2010, we extrapolate the magnetic field to the upper 
atmosphere from the photosphere. By calculating the time variation of magnetic field, we 
can get the induced electric field. The derived induced electric field can reach a value 
of 10 2 V/cm and the average electric field has a maximum point at the layer of 360 km 
above the photosphere. The Monte Carlo statistics method is used to compute the triple 
integration of induced electric field. 
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1 INTRODUCTION 

Electric field in solar atmosphere plays an important role in heating plasma, and accelerating and trans- 
porting charged particles (Priest & Forbes 2000). At the same time, its distribution provides rich in- 
formation of solar flare, as well as other dynamic solar activities. The simultaneous determination of 
electric and magnetic field vectors make the estimation of Poynting flux of electromagnetic energy en- 
tering the corona and the flux of relative magnetic helicity available. However, the determination of 
electric field is quite harder and less than that of magnetic field in the solar atmosphere. 

Although the solar activities are dominated by magnetic filed, and much progress was made in this 
aspect in the past decades, there are still many points out of our understanding, such as the physical 
mechanisms of flares and filament eruptions. It is time to figure out that electric field, its magnitude, 
geometry, time-dependence and especially spatial distribution do provide us powerful tools to probe in 
solar activities where particle accelerating and energy release are believed to occur. 

Wien (1916) is the first author who pointed out the possibility of measuring electric field of solar 
plasma and measured the motional electric field using Stark effect. Around 1980s, some attempts to 
measure the electric field with helium and silicon spectra had been made, which suggested an electric 
field of 700 V/cm (Davis 1977) and 300 V/cm (Jordan et al. 1980) respectively. Even so, Moran & 
Foukal (1991) pointed out that direct measurement of the electric field by Stark effect is hard to operate 
with the low sensitivity. 

Because there was not an efficient and reliable way to get the direct electric field of solar atmosphere, 
researchers thought out many indirect methods, which either explicitly or implicitly invoke the ideal 
MHD assumption, E = —v x B/c. Tracking method and inductive method are two classes of such 
techniques. 
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According to the ideal MHD equation above, the magnetic field is known from the vector magne- 
togram, so both of two classes have focused on determining the velocity vector. The tracking method, 
developed by November & Simon (1988), compute velocity through a cross-correlation function that 
depends on the shift of feature points between two images. Although tracking methods are simple and 
robust, they also suffer from some shortcomings, for example, this technique is actually two-dimensional 
without the vertical component. Inductive methods , first developed by Kusano et al. (2002), improves 
the result of tracking methods with a solution to the vertical component from the magnetic induction 
equation and derives a three-components velocity vector from a sequence of vector magnetograms. 
Since the seminal work of Kusano et al. (2002), several techniques have been developed to determine 
velocity from vector magnetograms, and Welsch et al. (2007) provided detailed tests and comparisons 
of these techniques. 

Poletto & Kopp (1986) derived the maximum electric field of 2 V/cm in a large two-ribbon flare 
using the reconnection theory of Priest & Forbes (2000) where they used a very simple relationship 
between the electric filed along the current sheet and the observable velocity and magnetic field. With 
the similar theory, Wang et al. (2003) discovered two stages of electric evolution in another two-ribbon 
flare: electric field remains 1 V/cm averaged over 20 minutes as the first stage, and are followed by 
0.1 V/cm in the next 2 hours. Qiu et al. (2002) worked on an impulsive flare with high cadence Ha 
observations at Big Bear Solar Observatory (BB SO) and estimated the maximum electric field to be 90 
V/cm. 

Recently, Fisher et al. (2010) pointed out a way to compute electric field from a sequence of the 
vector magnetograms by using the Faraday's law and showed that it is possible to derive electric field 
whose curl is the time derivatives of three-components of B. The main problem for these authors is the 
non-unique solution of Faraday's law. 

In this paper, our goal is to present a method to indirectly compute the three-components induced 
electric field distribution in solar atmosphere through evolution of vector magnetic field. We first extrap- 
olate the magnetic field from photosphere to corona, and then calculate the induced electric field from 
the time variation of magnetic field. 

In solar plasma, there are three kinds of electric field: static electric field, induced electric field and 
motional electric field. Here we only focus on induced electric field which is caused by the change of 
magnetic field. The static electric field which is triggered by existence of charged particles accumulation 
and the motion of plasma in the direction perpendicular to magnetic field are not considered due to the 
screening of plasma in the solar atmosphere. 

The paper is structured as follows. In section 2 we extrapolate vector magnetic field distribution 
from magnetogram; in section 3 we do the computing of electric field from the extrapolated magnetic 
field; and next in section 4 we give an examination by use of the observational vector magnetograms 
from HMI/SDO. Conclusion and discussion are given in section 5. 

2 EXTRAPOLATION OF MAGNETIC FIELD 

At present, although many attempts had been done to estimate the coronal magnetic field (House 1977; 
Arnaud & Newkirk 1987; Judge 1998; Judge et al. 2001), the reliable information about magnetic field 
is available only for the photospheric level. Just similar to the electric filed, we currently and in the near 
future have to face the fact that the direct measurements of magnetic field in the global solar atmosphere 
are still unavailable. 

There are many methods to extrapolate magnetic field from the photospheric magnetic field under 
the assumption that the magnetic field is nearly force-free. Force-free magnetic field of solar atmosphere 
must satisfy the following equations: 



jxB = 0, 



(1) 



V x B = aB, 



(2) 
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V • B = 0, (3) 

where a is a scalar function of position and time. The above equations imply nonexistence of Lorentz 
force and a being a constant along the magnetic field line. The equations represent a potential field if 
a = 0, a current-carrying linear force-free(LFF) field if a = constant, and a general nonlinear force-free 
field(NLFF) field if a = f(r) respectively. 

The extrapolations of potential and linear LFF field are maturely developed. Potential field and 
LFF field can be determined directly from the line-of-sight (LOS) component of magnetic field(e.g. 
MDI/SOHO) as input, and a has to be computed in LFF field from some additional data(Chiu & Hilton 
1977; Seehafer 1978; Alissandrakis 1981; Gary 1989). 

For the NLFF field, several methods have already been pointed out: Grad-Rubin method (Sakurai 
1981, MHD relaxation method (Chodura & Schlueter 1981 Roumeliotis 1996) and the optimization 
method (Wheatland et al. 2000). The last one will be used in this paper. 

In the optimization approach, Wheatland et al. (2000) defined a quantity L: 



L = I [B- 2 \(\7xB)xB\ 2 + \\/.B\ 2 }dV, 
Jv 



(4) 



where B is defined in a volume V. If L is decreased to zero, equations (l)-(3) are fulfilled, then the 
field is force-free in the volume V. In order to reduce L, B needs to be evolved like: 

dB 

1H=» F - (5) 

Wheatland et al. (2000) provided a testing to optimization method. Inhester & Wiegelmann (2006) 
provided detailed comparison of optimization method and Grad-Rubin method by implementing these 
two algorithms and comparing the performance. In addition, Liu et al. (201 lb) used two semi-analytical 
solutions of force-free fields to test two other NLFF extrapolation methods: boundary integral equation 
(BIE) method developed by Yan & Sakurai (2000) and approximate vertical integration (AVI) method 
developed by Song et al. (2006). 

Wiegelmann (2004) improved the optimization approach by showing how the magnetic field can be 
reconstructed only from the button boundary and developed a code which will be used later. 

Moreover, although the NLFF field model is widely-used (Rgnier 2004, Wiegelmann et al. 2005, 
Schrijver et al. 2008), a joint study of De Rosa et al. (2009) concluded that a successful application of 
NLFF field extrapolation should satisfy several requirements. Recently Wiegelmann et al. (2012) offered 
a detailed discussion on this problem and proved that their results fulfill these requirements. The results 
of nonlinear force-free modeling should be used with some caution. 

3 CALCULATING INDUCED ELECTRIC FIELD DISTRIBUTION 

In solar atmosphere, the static electric field is neglected due to the plasma screening. We consider the 
case that there is no accumulation of charged particles in the solar atmosphere, and the electric field is 
mainly generated from the time variation of magnetic field. That is, 

V • E = 0, (6) 

and 

V x E = fl (7) 

where 

dB 

at 

We could derive electric field directly from the time variation of magnetic field (Batchelor 2000): 
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In order to get E from equation (9), we introduce the Monte Carlo method which is a statistical 
simulation method to solve the triple integral problem. This method can be used to approximate the 
involved integral if the precise value of integrals is not important and its estimated value is enough, or 
if the precise value is unable to get. 

If f(x,y,z) is continuous function on domain D, and g(x,y,z) is a probability density function such 

that 



g(x,y, z)dxdydz = 1, (10) 

id 

(xi, yi, Zi)(i = 1,2, ...N) are a sequence of random numbers that fall in domain D, according to the 
theorem of Monte Carlo method, when N is large enough, we have 

lllj( X ,y, Z ) dxdydz ^^i^^. (11) 



If g(x,y,z) is constant, equation (10) becomes 

f 



1 



and 

D N 

f(x, y, z)dxdydz w — ^ f(xi,yi, Zi). (13) 
. By applying Monte Carlo method to equation (9), we get the final equation of E as 

E -J^l^ R3 ( 14 > 

i=l 



4 THE INDUCED ELECTRIC FIELD IN NOAA AR11072 



To implement the method we described in the previous sections, we provide an example here. At first 
we extrapolate magnetic field of photosphere from magnetograms, then we compute the time different 
of two magnetic field distributions ft = —dB/dt, at last we use Monte Carlo method to calculate the 
electric field distribution. 

In this example, we use observational data from HMI/SDO which provides high spatial and temporal 
resolution vector magnetograms(Schou et al. 2012). HMI provides continuous vector magnetograms at 
a 12-minute cadence and released in the past two years several data series of cutouts of the original 
full disk images. The full released data summary can be found at the SDO Joint Science Operations 
Center web page (http://jsoc.stanford.edu). From the available data releases, we chose the 6-day cutouts 
of 5 1 2 x 5 1 2 pixels for NOAA Active Region 1 1 072 from May 20th 20 1 to May 26th 20 1 0. 



Table 1 Maximum and average value of B at different layers 



z 

(km) 


max(\B x \) 


max(\B y \) 


max(\B z \) 
(G) 


mean(B x ) 


mean(B y 


) 


mean(B z ) 


360 


1 x 10 3 


8.4 x 10 2 


8.3 x 10 2 


4 x 10" 1 


-1.1 x 10" 


-l 


-3.2 


720 


1.1 x 10 3 


7.8 x 10 2 


1.1 x 10 3 


7.6 x 10" 1 


-1.2 x 10" 


-l 


-3.4 


1080 


7.2 x 10 2 


1.1 x 10 3 


2.7 x 10 3 


6.4 x 10" 1 


-9.6 x 10" 


-2 


-3.2 


1440 


8.2 x 10 2 


1 x 10 3 


4.5 x 10 2 


4.5 x 10" 1 


-3.3 x 10" 


-2 


-3.3 


1800 


1.9 x 10 2 


4.2 x 10 2 


8.6 x 10 2 


3.9 x 10" 1 


-4.4 x 10" 


-2 


-3 


2160 


2.7 x 10 2 


2.7 x 10 2 


8.1 x 10 


3.2 x 10" 1 


-3.3 x 10" 


-2 


-3 
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Fig. 1 Extrapolated results of a magnetogram for Active Region 11072 observed from 
HMI/SDO at 12:12 on May 20th, 2010. The size of every image is 10.8 x 10.8 Mm 2 , and 
the value of z means the distance from photosphere in the unit of km. Arrows show directions 
and amplitudes of B x and B y , while the background image shows the amplitude of B z 



We run the code developed by Wiegelmann (2004) (optimization method) to extrapolate magnetic 
field from vector magnetograms, and both potential and NLFF field are generated, Liu et al. (2011a) 
provided detailed comparison of NLFF Field and Potential Field. We extrapolate magnetic field to a 
volume of 512 x 512 x 9 pixels and the distance between two adjacent pixels is 0.5 arc seconds. Figure 
1 shows the extrapolated NLFF field of one magnetogram. 

The horizontal cuts of the lower six layers consisting 30 x 30 pixels are shown in Figure 1 . The value 
of z indicates height from the photosphere in unit of kilometer, arrows show directions and amplitudes 
of B x and B y , while the background image shows the amplitude of B z . 

Table 1 shows the maximum and average value of the three components of the extrapolated magnetic 
field, and the first image of Figure 5 shows the average of absolute value of magnetic field at six layers. 
We show in Figure 2 the difference between two extrapolated results, that is, B2 - B\ where £?i is the 
magnetic field in Figure 1 . 



Table 2 Maximum and average value of A£? at different layers 



z 

(km) 


max(\AB x \) 


max(\AB y \) 


max(\AB z \) 
(G) 


mean(AB x ) 


mean(ABy) 


mean(AB z ) 


360 


1.2 x 10 3 


1.4 x 10 3 


1.2 x 10 3 


1.5 x 10" 1 


1.7 x 10" 1 


4.3 x 10" 


-2 


720 


1.1 x 10 3 


1 x 10 3 


1.3 x 10 3 


2.2 x 10" 1 


2.7 x 10" 3 


4.1 x 10" 


-2 


1080 


7.4 x 10 2 


1.1 x 10 3 


2.7 x 10 3 


-1.7 x 10" 1 


1.6 x 10" 1 


4.2 x 10" 


-2 


1440 


8.1 x 10 2 


1 x 10 3 


4.8 x 10 2 


-8.9 x 10" 2 


-6.4 x 10" 3 


3.9 x 10" 


-2 


1800 


1.9 x 10 2 


4.2 x 10 2 


8.7 x 10 2 


-8.8 x 10" 2 


6.8 x 10~ 2 


3.7 x 10" 


-2 


2160 


2.7 x 10 2 


2.6 x 10 2 


8.9 x 10 


-9 x 10~ 2 


3.7 x 10~ 2 


3.6 x 10" 


-2 
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Fig. 2 Time variation of magnetic field extrapolated from two continuous vector magne- 
tograms, the first one is observed at 12:12 on May 20th, 2010 by HMI/SDO, and the other is 
observed 12 minutes later. 



Table 2 shows the maximum and average value of the difference of the extrapolated magnetic field, 
and the second image of Figure 5 shows the average of absolute value at six layers. Monte Carlo method 
requires that ft = — ^ in equation (9) and (14) is spatial continuous function, but our observed and 
extrapolated Q is discrete with spatial distance of 0.5 arc seconds. Thus, we linearize Q in unit of 0.5 
arc seconds to reconstruct a continuous function. 




Fig. 3 The ft at the point A(x,y,z) is linearized from 8 points next to it. The position of these 
adjacent points is (xi, yi, Zi, i = 1...8), and every edge of the cube represents 0.5 arc seconds. 
A has two projections in the upper and lower surfaces where the Q should be computed at 
first. 
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We use the linear algorithm below to compute Q(x, y, z u p) of the upper projection in Figure 3: 

Vback, Zup) = to(xieft, Vback, Z up ) X {x right x)/0.5) + ^(x right: Vback, Z u p) X (x - Xi e f t )/0.5, 
Q(X, y forth, Z up ) = Q(xi e ft, y forth, Z up ) X (x righ t x)/0.5) + Q(x righ t , 2/ forth, Z up ) X (x - Xi e f t )/0.5, 

Q(x, y, z up ) = Sl(x, yback, z up ) x (y fort h - 2/)/0.5 + Sl(x, y forth, z up ) x (y - y bac k)/0.5, 

where 

Xleft = Xi = X4 = X 5 = X 8 , X r ight = X 2 = X 3 = X 6 = X 7 , 

Vback =2/1=2/2 = 2/5= 2/6, V (forth) =2/3 = 2/4 = 2/7 = 2/8, 

Zdown — Z§ = Zq = Zj = Zs, Z up = Z\ = Z 2 = Z3 = Z4. 



z=360km 


"Dili. ♦ 1 r^y 1 r ^ ; " 
:::;;^/S::tt:::: 






. f . , , , , , „^ 

* * * ^ ^ ':*> ," 

:::: 


/ 1 w - 


t * t f ^ , v * 

- * * * * * * ~.g# / * * * ' ' " 






- - a. 

z=1 800km 


- ' ' ' * - i J, * * " * ' ' * * - 


z=2160km . ; " 



Fig. 4 Computed "vector electrograms" (10.8 x 10.8 Mm 2 ) of different depth. Arrows show 
directions and amplitudes of and E y , while the background image shows the amplitude of 

E z 



Similarly, we can get Q(x,y, Zdown), together with Q(x,y, z up ), and we compute the linearized 
Q(x, y, z) from 

Sl(x, 2/, z) = Q(x, 2/, z down ) x (z up - z)/0.5 + Sl(x, y, z up ) x [z - z down )/0.b. 

For saving time, we only generate 10,000 random points in the cube of 512 x 512 x 7, and calculate 
the distribution of induced electric field depending on the equation (14). One of our result are showed 
in Figure 4. 

Table 3 shows the maximum and average value of the calculated three components of the electric 
field and the third image of Figure 5 shows the average of absolute value of electric field at six layers. The 
average of absolute electric field reaches maximum point at the layer of 360 km above the photosphere. 
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Table 3 Maximum and average value of E at different layers 



z 


max(|^|) 


max( \E y \) 


max(|^|) 


mean( J E cc ) 


mean(i^) 


mean(i^ ) 


(km) 


(V/cm) 






360 


5.5 x 10 2 


3.9 x 10 2 


2.5 x 10 2 


3.5 x 10~ 3 


1.6 x 10" 3 


1.1 x 10" 3 


720 


3.6 x 10 2 


1.7 x 10 2 


4.4 x 10 2 


9 x 10~ 4 


2.3 x 10~ 4 


-1.4 x 10" 3 


1080 


4.6 x 10 2 


2.3 x 10 2 


1.6 x 10 2 


2.1 x 10" 3 


1.2 x 10" 3 


7.9 x 10" 4 


1440 


8.1 x 10 


2 x 10 2 


1.1 x 10 2 


-2.1 x 10" 4 


-7.5 x 10" 5 


-3.5 x 10" 5 


1800 


1.2 x 10 


1.1 x 10 


1.8 x 10 


-2.7 x 10~ 4 


-4.3 x 10~ 4 


-1.6 x 10~ 4 


2160 


2.2 x 10 


1.2 x 10 


9.3 


-2.2 x 10" 5 


-4.1 x 10" 4 


-1.2 x 10" 4 



5 SUMMARY AND DISCUSSION 

In this paper, we describe and implement a new method to calculate the distribution of induced electric 
field in solar atmosphere using a sequence of vector magnetograms as an input. 

We first introduce several extrapolation methods of magnetic field and do a simple comparison of 
these methods, then we choose the optimization method in our example to extrapolate the magnetograms 
observed by HMI/SDO from photosphere to corona. After that we derived a special solution of electric 
field in the form of triple integral. 

To solve the triple integral problem, we utilize the Monte Carlo statistical method to get a new 
equation of electric field. As this method require the continuous function, we linearize ft which is 
originally spatial discrete function. The similar linearization has been used to compute time variation 
of magnetic field from the magnetograms of 12-minute time resolution, that is, we assume that B go 
through linear change in this 12-minute interval. 

Through the derivation, it is proved that as long as boundary condition (Batchelor 2000) is fulfilled, 
we can obtain three-component electric field of solar atmosphere only from vector magnetograms. In our 
example for NOAA AR1 1072, the result shows that intensity distribution of induced electric field varies 
at different layer: it reaches a value of 10 2 V/cm and the average electric field has a maximum point at the 
layer of 360 km above the photosphere. 1 However, there are several shortcomings of this method needed 
to be figured out. At first, the boundary condition (Batchelor 2000) is not strictly satisfied, because in 
the quiet area, time variation of magnetic field is small but not zero. Secondly, the temporal and spatial 
resolution of the HMI/SDO vector magnetogram is still not high enough to provide continiuous time 
and spatial sequences, linearlization methods have to be used twice in our calculation: one is to compute 
dB I dt which assume magnetic field linearly change at this 720s interval, and the other is to satisfy 
the requirement of Monte Carlo method which needs spatial continuous function as input. Thus our 
calculation process might not fully reflect the real situation of the sun. Finally, we use Monte Carlo 
method to compute the triple integration of a large volume which is a time-consuming job, and we have 
to adapt a small sample of only 10,000 random numbers which is far away from large enough. 
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Fig. 5 The average of absolute intensity of B, AB and E at six layers. 



